Controlling open quantum systems using fast transitions 
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Unitary control and decoherence appear to be irreconcilable in quantum mechanics. When a 
quantum system interacts with an environment, control strategies usually fail due to decoherence. 
In this article we implement a time-optimal unitary control protocol suitable for quantum open 
systems. The method is based on succesive diabatic and sudden switch transitions in the avoided 
crossings of the energy spectra of closed systems. We show that the speed of this control protocol 
meets the fundamental bounds imposed by the quantum speed limit, thus making this scheme ideal 
for application where decoherence needs to be avoided. We show that this method can achieve 
complex control strategies with high accuracy in quantum open systems. 
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Introduction - Quantum control is a fundamental 
goal in different areas, including physical chemistry, 
nanoscience and quantum information processing pj[3]. 
In fact, as the manipulation of quantum systems will be 
the basis for future technological applications, the deve- 
lopment of different control strategies is a task of major 
interest nowadays. 

An unitary control scheme consists of engineering the 
time dependence of one or several external parameters 
(i.e. electric fields) which manipulate the temporal evolu- 
tion of a quantum system in order to obtain a desired tar- 
get state. But, real systems interact with their environ- 
ment to a greater or lesser extent. No matter how weak 
the coupling between the system and the environment is, 
the evolution of an open system will be characterised by 
nonunitary effects like decoherence and dissipation [4]. 
Decoherence, in particular, is a quantum effect whereby 
the system loses its ability to exhibit coherent behavior. 
It is a serious obstacle in quantum information process- 
ing in general, and in quantum control in particular [5]. 
Control faces major difficulties when dealing with open 
systems, as the control field cannot fully compensate the 
dissipative and decohering effects from the environment. 
Facing this threat requires a control protocol not only to 
assure high fidelity but also to be implemented in a time 
span shorter than the decoherence time [6j [7] . Neverthe- 
less, quantum mechanics imposes fundamental limits on 
the maximum rapidity with which a state can evolve to 
an orthogonal state. Known as the quantum speed limit 
(QSL), it is a physical bound that any control strategy 
should take into account [8HT0]. 

In this paper we implement an efficient method for 
controlling a quantum system in contact with an envi- 
ronment. The method is based on the knowledge of the 
spectrum of the closed system as a function of a suitable 
external parameter and requires that the system behaves 
locally -near avoided crossings- as the Landau- Zener (LZ) 
two level model JTTJ [12]. Although this characteristic 
may seem rather restrictive, it is, in fact, a general prop- 



erty of systems with interaction between its energy levels 
p~3| H4] , at least in the low energy region. We apply a 
well-defined series of fast (diabatic) and sudden (step- 
like) variations of the control parameter, which allows us 
to travel through the state space of the system and reach 
the desired target state. The speed of this control proto- 
col meets the fundamental bounds imposed by the QSL, 
thus making this scheme ideal for application to open 
systems. We show that this method is very successful to 
accomplish ambitious control goals in systems exposed 
to an important environmental decoherence. Addition- 
ally, this method avoids the need of using optimization 
algorithms in order to find a time-optimal control func- 
tion, and thus various potentially encumbered subjects 
such as monotonicity, convergence and performance of 
the algorithm. 

Control strategy and QSL - The building block of this 
strategy is the LZ model for a quantum system with two 
interacting levels. In the diabatic basis {|0) , |1)} the 
hamiltonian of this system can be written 

where A is a constant and we set ft = 1 along the Letter. 
The eigenvalues of Hlz form an hyperbola in the (A, E) 
plane (as shown in Fig. 1 (a)), whose vertex represents 
an avoided crossing (AC) with an energy gap A. The 
eigenvectors of the hamiltonian form the so-called adia- 
batic basis, {|</>o) , where the notation makes ex- 

plicit the asymptotic correspondence between both basis 
when A — ^ — oo. When A — > -f oo, this correspondence 
is exchanged. The classic LZ theory describes the tran- 
sition probability of the system when the initial state is 
oo)) = |0) and the parameter A is sweept linearly in 
time, i.e. X(t) = v t, yielding the famous LZ formula 

P l( ^oo) = l-exp(-^). (2) 
This result defines a critical velocity v c = ^h- which de- 



2 



termines two limiting control scenarios. In the first place, 
we have the fast diabatic (D) transitions, in which v ^> v c 
in such a way that P\{t — »• oo) c± 0, leaving the initial 
state unchanged. On the other hand, the adiabatic (A) 
transitions, in which v <C v c and thus P\{t — » oo) ~ 1, 
take place when the state evolves slowly following the 
adiabatic curve and finishes in state |1). For a closed 
quantum system in which the transitions between neigh- 
bouring levels are well described by the LZ model, this 
scheme provides a quantitative binary (A-D) recipe for 
determining the appropiate velocity of the driving field 
at each of the avoided crossings (ACs) jT5HT7j . However, 
in the presence of an environment, decoherence is bound 
to act within the long periods of time required by adia- 
batic transitions, rendering the effective dynamics of the 
system non-unitary and thus preventing the system from 
being controlled. 

To overcome this problem we propose an alternative 
method in which we use the so-called sudden-switch tran- 
sitions p~6j [18]. Consider the LZ system prepared ini- 
tially in the state ^(O)) = |0). We now consider \{t) 
to be a piecewise constant function with initial value 
A(0) = — Ao, with Ao ^> A/|a| [13 . In this way, the ini- 
tial state is approximately an instantaneous eigenstate of 
the Hamiltonian Hlz(0). If now X(t) undertakes a sud- 
den variation to A = and the system is left to evolve 
for a time T, the final state is given by 

|^(r)) = cos(|r)|o> + sin^r)|i>. (3) 

It is clear that choosing T = f yields |^(T)> = |1). The 
final step in the evolution is a second sudden switch 
of the control parameter from to +Ao. Once again, 
the instantaneous eigenstates of the hamiltonian will be 
approximately those of a Z: in such a way that \ijj(t)) will 
be a stationary state for t > T, since \^(T)} = |1). In 
sum, we have driven the system from |0) to |1) in a time 
T = J with a probability of 1. In this way, this scheme 
represents a sort of shortcut to adiabaticity, as the final 
state of the system is the same as if the parameter \(t) 
had been modified adiabatically. The total evolution 
time has been dramatically shortened, since an adiabatic 
passage through the AC requires a total time much 
larger than a critical time of the order of A -2 (recall the 
critical velocity expression from the LZ formula). More- 
over, the total evolution time T is precisely the shortest 
possible time in which this two level system can change 
its state, as a consequence of the time-energy uncertainty 
principle. This can be seen as follows. Bhattacharyya 
[9] derived from the Mandelstam-Tamm inequality [19] 
the expression r > arccos(v^r)/AiI, in which P T is 
the quantum non-decay probability, i.e., P T = | ((/>|(/> r ) | 2 , 
where \(j> t ) = exp (— ij-i) In the context of control 
theory, r can be chosen as the total time of the control 
protocol (for the two level system), in such a way that 
\[P r = \(ip(0)\i/jgoai)\- Evaluating this result for the 



hamiltonian ([I]), with |^(0)) = |0) and \i/j g0 al) — |1) we 
get the minimum possible r to be tqsl = tt/2AHq. 
AHo can be calculated straightforwardly recalling from 
^ that H = H(X = 0) = f a x . This yields the result 
tqsl = J = T, i.e., exactly the total evolution time of 
the sudden-switch transition, which makes it the best 
possible candidate for replacing slow adiabatic driving 
at an avoided crossing, as its time-optimality favors the 
prospect of avoiding decoherence. 




Figure 1. Fidelity as function of the dimensionless time param- 
eter r = tA for sudden-switch variations of the control parameter, 
(a) Energy spectrum as a function of parameter A for the LZ hamil- 
tonian The asymptotic correspondence between the diabatic 
and adiabatic states has been made explicit, (b) Fidelity as a func- 
tion of r for adiabatic passages through the AC. In both fidelity 
plots, the solid line corresponds to the closed-system case, and the 
dashed curves to the system coupled to a bosonic environment. 
The temperature of the bath was fixed at T = 10A. 



Controlling the dissipative LZ model - We now show 
how our control method gives satisfactory results when 
applied to the LZ hamiltonian coupled to a bosonic en- 
vironment. We propose a master equation approach in 
the weak coupling and high temperature limit, which be- 
comes independent of the particular election of time de- 
pendence of the control parameter. Assuming bilinear 
coupling of the form a z (g) ^ZjCjqj (the g/s being the 
position operators of the environmental oscillators), the 
master equation for the reduced density matrix p(t) re- 
sults @] 



p{t) = -i(H\t)p(t)-p{t){H f )\t)) 



7oT 



(4) 



where 70 is the (dimensionless) coupling constant be- 
tween the system and the environment at equilibrium 
temperature T (we set the Boltzman constant ks = !)■ 



H'(t) = X(t)a z + f (1 



a T is the renormalized non- 
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hermitian effective hamiltonian of the system. This dif- 
ferential equation describes the non-unitary evolution of 
the reduced density matrix p(t) of the system which, 
in the small coupling limit, will be mainly due to the 
double-commutator term. We wish to study how this 
non-unitary dynamic affects the fidelity of the control 
scheme, defined in this context as F(t) = Tr (p(t)p goa i). 

It has been shown [20j [21] by different analytical ap- 
proaches that while diabatic (p g0 ai = |0)(0|) transitions 
in a dissipative LZ do not seem to suffer decoherence 
and thus achieve high fidelity, adiabatic (p goa i — 
variations of the control parameter render poor final 
state fidelity even for weak coupling. This can be seen 
in Fig. [I] inset (b), where the fidelity as a function of 
time is plotted for different values of 70. As A(t), being 
sweept slowly, reaches the AC, the system becomes 
extremely sensitive to decoherence and rapidly becomes 
mixed. The fidelity then fails to achieve the desired 
value of 1. For sufficiently large 70, the state losses 
its purity and evolves to p = I, for which J-fi na i = §• 
However, the desired transition to p goa i — is 
indeed achieved with high probability if the system is 
driven by the sudden switch method. In Fig. [I] we show 
how the fidelity in this case evolves favorably even in the 
presence of the environment, making this method much 
more robust under the action of external influences, 
due to its time-optimality. Moreover, the difference 
in the total evolution time between both methods is 
remarkable, as the adiabatic driving takes up to three 
orders of magnitude more time than the sudden-switch 
method. In sum, this represents a time-optimal, rela- 
tively decoherence-resistant control recipe for achieving 
"adiabatic" (Pi ~ 1) transitions at an AC in the energy 
spectrum of a system. Together with the diabatic 
transitions, this updated binary control protocol (which 
we will refer to as S-D) renders a shorter control time 
and high fidelity at each AC, even in the presence of a 
dissipative environment. 

Controlling a dissipative multilevel system - We will 
show now how to apply our control method to a more 
complex multilevel system. Consider two interacting 
spin-^ particles £4 and <Sg, with dynamics given by the 
hamiltonian 

H = A A a M + cos *i B) + A B ff W + 5 , (5 ) 

where cr\ represents the i-th Pauli operator in Hilbert 
space associated to particle X and all energy parameters 
are fixed but uj% = A A ^4. In Fig. [2] we show the energy 
spectrum associated to hamiltonian ([5| as a function of 
A. This spectrum resemble the one corresponding to re- 
alizable systems like a qubit in a cavity or a Josephson 
junction qubit [22]. The spectrum is found to have four 
ACs where the levels interact; two of them are located 
at A = while the other two are placed symmetrically at 
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Figure 2. Spectrum of the four-level system |5f, i. e., di- 
mensionless energy e — | as a function of parameter A. The 
remaining parameters of the hamiltonian were set as A^ = 50 and 
A B = 0.05A A . The side ACs have a gap of A = 2.43, while the 
central ACs have wider gaps of size 2A. The capital letters refer to 
the diabatic state corresponding to each energy branch at A — y ±00. 
The symbols ■ and refer to the initial and target states, respec- 
tively, which are joined by a path following the energy levels. Full 
segments of the path symbolize diabatic variations of the control 
parameter, while dashed lines refer to sudden-switch transitions. 



A = ±A C . We choose the parameter S = 0.5A^4, in such a 
way that two of the ACs have a gap A and the remaining 
two have a wider gap of 2 A. The correspondence between 
the instantaneous eigenstates of the hamiltonian H(X) 
and the diabatic basis {\4>a) ? \<I>b) •> |0c) •> I&d)} has been 
made explicit in the figure. Since we are specially inter- 
ested in studying how the fidelity of the method evolves 
in the presence of decoherence, we once again consider a 
thermal bosonic environment coupled to the subsystem 
bilinear ly as (g) ^ZjCjqj. The master equation 
for the reduced density matrix of the composite system, 
in the small coupling and high temperature limit, turns 
to be identical in form to Q, where a z and A are to be 
replaced with and A^. 

We now pose the problem of starting from an initial 
state i/)(0) = 4>a and driving the system to ip goa i — 0c- 
For that purpose, we design a path in the energy spec- 
trum, which can be seen in Fig. |2) This path requires the 
control parameter to be sweept following the sequence D- 
S-D-S (where D stands for a diabatic transition and S for 
sudden-switch transition) in order to achieve the desired 
navigation of the spectrum. In Fig. [3] we show the time 
dependence chosen for A(t) in order to achieve this evo- 
lution. The high-slope linear segments correspond to di- 
abatic transitions, and two of the the constant segments 
have a lenght tt/A and 7r/2A, respectively, as required 
by the sudden-switch method. 

In Fig. [3](b) we show the evolution of the probabilities 
of finding the system in each of the various diabatic 
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Figure 3. (a) Control parameter A as a function of the dimension- 
less time parameter r = tA, designed to achieve the path displayed 
in the spectrum of Fig. [2] (see text for details), (b) Probabilities of 
finding the four-level system in each of the four diabatic states as 
a function of r. The dotted gray vertical lines indicate the times 
at which the two steps of the sudden-switches are performed, (c) 
Same as (b) but when the method is applied to the system coupled 
to an environment, with 70 = 10 -3 . The final fidelity is ~ 0.85. 
The dashed curve represents the evolution of the purity of the state. 
The parameters of the system are the same as those specified in Fig. 
[2] The temperature was set to T = 20 A. 



parameter is sweept linearly with high speed, after which 
the state of the system is still Finally, a second 

sudden switch evolves the system from (/>£> to <pc, now 
in a time 7r/2A. Each step of the protocol is successful, 
yielding an elevated final fidelity of F(T) c± 0.99. More- 
over, the total time of evolution has barely added up to 
27r/A, which represents an excellent time performance. 
The results of applying the same control protocol to the 
system coupled to the environment (70 7^ 0) are shown 
m Fig. [3] (c). In this case, the non-unitary fashion of 
the evolution manifests itself through the decay of the 
purity of the state in time, although the rapidness of 
the control protocol allows it to act in a time interval 
shorter than the decoherence time, making the system 
follow the desired path and consequently achieving 
a fairly high fidelity. Moreover, as most of the total 
evolution time is due to the sudden-switch transitions, 
which are time-optimal, its bound to expect that there 
is no room for improving the final state fidelity with 
this or any other kind of unitary control. In Fig. [4] 
we show how the final fidelity T varies as the coupling 
with the environment is increased. There it can be 
appreciated how the S-D method yields a fidelity of 
over 0.8 for 70 < 10 -3 , and becomes unsuccessful for 
sufficiently large 70, as can be expected due to the 
presence of an strongly decohering environment. For 
comparative purposes, a similar plot is shown with the 
results obtained applying the previous A-D method to 
the four-level system, where we can see that the old 
method fails even for small coupling. 
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Figure 4. Final fidelity F(T) as a function of the dimensionless 
coupling constant 70, for the four-level system, using S-D and A-D 
control methods. 



states, when the driving field of Fig. [3] (a) is applied to 
the closed system (70 = 0). The fidelity in this case is 
F{t) = Pc(t), i- e -> the black curve in the figure. There, 
it can be seen how the state remains unchanged, as 
expected, after the first diabatic passage. Then, as the 
control parameter undertakes the first sudden-switch, 
the state evolves to <\>b in a time tt/A. Next, two 
sucesive diabatic transitions take place when the control 



Final Remarks - We have implemented an efficient 
method to control the state of a quantum open system. 
The method is based on the navigation in the energy 
spectrum of the closed system, using fast variations of 
a control parameter (which are already familiar schemes 
used in NMR, quantum optics, and cavity and circuit 
QED, for engineering quantum dynamics). The success 
of this method to overcome the adverse influence of de- 
coherence relies in the speed in which the transitions are 
performed. We show that the sudden transitions occur 
in the quantum speed limit. The method was succesfully 
applied in a system with four energy levels to reach de- 
sired target states that lie far in the spectrum from the 
initial state. We stress that our proposed control function 
has a simple analytical form even for complex multilevel 
systems, as opposed to those obtained through optimiza- 
tion methods (see for example [23]). This fact makes our 
method very attractive for future experimental applica- 
tions where several of the requirements of the method are 
fulfilled, but decoherence appears as an obstacle for the 
quantum control. 
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